Subharmonic gap structure in short ballistic graphene junctions 
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We present a theoretical analysis of the current-voltage characteristics of a ballistic 
superconductor-normal-superconductor (SNS) junction, in which a strip of graphene is coupled to 
two superconducting electrodes. We focus in the short-junction regime, where the length of the strip 
is much smaller than superconducting coherence length. We show that the differential conductance 
exhibits a very rich subharmonic gap structure which can be modulated by means of a gate voltage. 
On approaching the Dirac point the conductance normalized by the normal-state conductance is 
identical to that of a short diffusive SNS junction. 

PACS numbers: 74.45.+c,74.50.+r,73.23.Ad,74.78.Na 



Introduction. - Recently it has become possible to ma- 
nipulate graphene, a single atomic layer of carbon, and to 
establish electrical contacts with rfci»2»&. This has opened 
the possibility of studying its transport properties, which 
are determined by the dynamics of two-dimensional Dirac 
fermions (massless excitations governed by a relativistic 
wave equation) 2 -'". Moreover, the fabrication of a new 
class of hybrid structures in which superconductors are 
coupled via graphene is now feasible^. It has already 
been shown theoretically that the nature of the low-lying 
graphene excitations might lead to unexpected features 
in transport properties of these heterostructures^. Thus 
for instance, Titov and Beenakker— have recently demon- 
strated that a SNS junction with an undoped strip of 
graphene as a normal region can sustain a supercurrent 
and it exhibits an unusual "quasi-diffusive" scaling of the 
dc Josephson effect. 

The supercurrent in graphene hybrid systems may 
nicely reveal the interplay of superconductivity and the 
relativistic quantum dynamics of the electrons in these 
carbon structures. However, a direct comparison of the 
supercurrent between theory and experiment is usually 
not straightforward. This is due to the fact that the 
superconducting phase is prone to both quantum and 
thermal fluctuations, which depend on the electromag- 
netic environment in which the junction is embedded 7 . 
Thus, unless the environment is carefully designed, which 
is not an easy task, the supercurrent might be greatly re- 
duced as compared with the theoretical predictions for 
idealized situations. For this reason, we propose to look 
at the current-voltage characteristics (I-V) of graphene 
SNS contacts, which do not suffer from the problems just 
mentioned. The main feature of the I-V curves of SNS 
junctions is the appearance of the so-called subharmonic 
gap structure, which consists of a series of conductance 
maxima at voltages 2A/ne, where n is an integer and 
A is the energy gap of the electrodes. This structure 
originates from the ocurrence of multiple Andreev reflec- 
tions (MARs) 8 . The microscopic theory of these tunnel- 
ing processes developed in the has been shown 



to describe accurately the I-V characteristics of super- 
conducting atomic point-contactsii. 

In this Communication we present a theoretical study 
of the I-V characteristics of SNS junctions, where the nor- 
mal region is a ballistic strip of graphene. We consider 
the experimentally relevant case of a short junction, in 
which the length of the strip is smaller than the supercon- 
ducting coherence length f. We show that the conduc- 
tance exhibits a pronounced subharmonic gap structure 
that can be tuned with a gate voltage. In the case of zero 
gate voltage, i.e. at the Dirac point, the I-V curves are, 
after rescaling, identical to those of a short diffusive SNS 
junction, which is a new consequence of the relativistic 
dynamics of the electrons in graphene. 

Model for the graphene junction.- We consider the sys- 
tem shown schematically in Fig. ^ (upper panel), where 
a graphene strip is coupled to two superconducting elec- 
trodes. We also assume that an additional gate electrode 
allows to control the carrier concentration in the normal 
part of the junction. We focus here in the case where the 
width of the graphene strip W is much larger than the 
junction length L (in this case the details of the micro- 
scopic description of the strip edges become irrelevant). 
Our main goal is the analysis of the superconducting 
I-V characteristics in the experimentally most relevant 
short-junction regime where L is small relative to £. In 
terms of energy scales, this condition requires A <C hv/L, 
where v is the energy-independent electron velocity in 
graphene. It has been shown that, as long as the normal 
transmission coefficients do not depend on energy on the 
scale of A, the transport properties of a short SNS junc- 
tion can be expressed as a sum of independent contri- 
butions from individual conduction modes 12 . Therefore, 
the transport properties depend only on the distribution 
of transmission probabilities Tj of these modes in the nor- 
mal state. For describing the transmission distribution 
of a graphene junction one can adopt the model intro- 
duced in Ref. [lj, which is briefly discussed in the next 
paragraphs. 

The authors of Ref. 13 considered the junction depicted 
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FIG. 1: (Color online) Upper panel: Schematic of a strip of 
graphene of width W , contacted by two superconducting elec- 
trodes (black rectangles) at a distance L. A voltage source 
drives a current through the strip. A separate gate electrode 
(not shown) allows the carrier concentration in the strip to 
be tuned around the neutrality point. Lower panel: Normal 
transmission of the conduction channels of the graphene junc- 
tion as a function of the transversal momentum in the limit 
W/L S> 1. The different curves correspond to different values 
of the dimensionless gate voltage k = e\V ga ,t a \L/hv. 



in Fig. ^ (upper panel). They determined the electron 
wave functions by solving the Dirac equation for mass- 
less fermions. In particular, they showed that assuming 
"infinite mass" boundary conditions at the strip edges 
y = and y = W, the transversal momenta are quan- 
tized as 

Qn = W n ^ + ' n = 0,1,2,..., (!) 

with n labeling the modes. Each mode has a twofold 
valley degeneracy. 

In the model of Ref. Il3l the gate voltage enters in the 
Dirac equation as an electrostatic potential. This poten- 
tial V(x) = V ga tc for < x < L determines the concen- 
tration of the carriers in the strip. The value V ga ± c = 
corresponds to charge neutrality, being the point where 
electron and hole excitations are degenerate (known as 
the Dirac point). The electrodes are modelled by tak- 
ing a large value V(x) — in the leads x < and 
x > L. (The parameter Voc drops out of the results, if 



\Voo\ ^> legate | •) The number of propagating modes N 
is given by N — Int (k^lY/ir + h), with el^l = hvk^. 
We are interested in the limit |14o| -> oo of an infinite 
number of propagating modes in the leads. By matching 
the solutions of the Dirac equation in the three regions of 
the junction one finds that the transmission probabilities 
at the Fermi level are given bjslii 

25 2 - 2(q n - k n f 2 

e k " L (q n - k n + i5) 2 + e- k ™ L (q n - k„ - iS) 2 

(2) 

with S = e\V gatc \/hv and k n = \/q 2 — S 2 for q„ > 5 
or k n — i\j ' 8 2 — q 2 for q n < 5. These transmssion co- 
efficients are plotted in Fig. ^ (lower panel) for different 
values of the dimensionless gate voltage k = e\V ga t e \L/Hv. 

The normal-state conductance Gn of the junction and 
the corresponding resistance i?N are given by 




n=0 



Superconducting I-V characteristics.- We now turn to 
the analysis of the current-voltage characteristics in the 
case in which the electrodes are in the superconducting 
state. As explained above, the current in the short- 
junction limit (L <C £) can be expressed as a sum of 
independent channel contributions as follows 

oo 

I(V,t) = 2j2HV,t,qn), (4) 

n=0 

where t is the time and I(V,t,q n ) is the single-channel 
current of a superconducting point-contact of transmis- 
sion T n {q n ) (cf. Eq. We compute this single-channel 
current using the Hamiltonian approach described in 
Ref. [Tot Let us remind that in the case of a supercon- 
ducting junction at finite voltage, the current oscillates 
on time as I(V, t) = J2 m Im(V) exp[im(j)(t)}, where 4>{t) 
is the time-dependent superconducting phase difference 
given by the Josephson relation dcf)(t)/dt = 2eV/h. We 
only consider here the dissipative dc current, which we 
shall simply denote as / from now on. Moreover, in the 
limit L <§; W considered here, one can replace the sum 
over n in the previous equation by an integration. 

We show in Fig. [5] the zero-temperature I-V charac- 
teristics for different values of the gate voltage. Notice 
that in the upper panel the current is expressed in an 
absolute scale, while in the lower one it is normalized 
by the normal-state resistance that scales as i?N oc L. 
At zero gate voltage, i.e. at the Dirac point, the cur- 
rent is formally identical to that of a short diffusive 
SNS junction^. This is easy to understand with the 
help of Eq. [5] From this expression one can show that 
for V^ a to = the transmission distribution adopts the 
form p(r) = {W/t:L)1/2t\/1 — r, which corresponds to 
the well-known bimodal distribution for diffusive wiresii 
Away from the Dirac point, the non-linearities of the I-V 
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FIG. 2: (Color online) (a) Zero-temperature current-voltage 
characteristics of a superconducting ballistic graphene junc- 
tion (length L short compared to the width W and the su- 
perconducting coherence length £), for different values of the 
gate voltage in the normal region (k = e\V ga ,t a \L/hv). (b) 
The same as in panel (a), but the current is normalized by 
the resistance of the junction in the normal state R^. 



curves are modulated by the gate voltage. In particu- 
lar, for gate voltages in which most of the open channels 
have a high transmission (see Fig. GJ>), the I-V curves 
are rather smooth (see for instance the curve for k = 4). 
Notice also that contrary to the normal state, the cur- 
rent at certain voltages below the gap is not neccesarily 
a monotonously increasing function of the gate voltage. 
This peculiarity is due to the fact that the subgap cur- 
rent is a very non-linear function of the transmission co- 
efficients. 

The non-linearities in the current are better observed 
in the differential conductance, which is shown in Fig. [3] 
The conductance exhibits a very pronounced subgap 
structure, although the maxima not always appear ex- 
actly at the voltages 2A/ne (n integer). As explained in 
the introduction, these peaks are due to the opening of 
new MAR processes at those voltages 8 '. The pronounced 
peaks (especially at 2A/e and A/e) come from the con- 
tribution of the low-transmitting conduction channels. 
This explains why this structure is more pronounced close 




FIG. 3: (Color online) Differential conductance correspond- 
ing to the I-V curves of Fig.|2] In panel (b) the conductance is 
normalized by the normal-state conductance Gn- As a guide 
for the eyes, the dotted vertical lines indicate the position of 
the voltages 2A/ne (n = I, 2, 7). 



to the Dirac point (see Fig. At very low bias, the 

conductance exhibits a square-root singularity (1/VV), 
which originates from the contribution of the highly- 
transmitting channels 12 . In an actual experiment, this 
singularity would be masked by the transition to the su- 
percurrent branch at low bias. 

Another interesting quantity for a possible comparison 
with experiments is the so-called excess current, I exc . 
At voltages eV ^> 2A the current can be expressed as 
I(V) — In(V) + Iexc, where If^((V) is the current in the 
normal state and I exc is the excess current, which is in- 
dependent of the voltage. In the case of a single-channel 
point-contact, the zero-temperature excess current as a 
function transmission coefficient t can be written asifi 



2eA 



h 1 - t 



1 



: In 



2(2 -t)VT=T \1-VT^P 



l + VT^r 



(5) 

Averaging this expression with the transmission function 
for the graphene junction, as indicated in Eq. one ob- 
tains the excess current shown in Fig.^J The gate modu- 
lation of I exc is very similar to the modulation of the crit- 
ical current found in Ref. 6, Again, at the Dirac point one 
recovers the result of a short diffusive SNS junctioni 2 ^, 
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eV L/hv 

gate 



FIG. 4: (a) Zero-temperature excess current as function of 
the gate voltage for a short ballistic graphene SNS junction, 
(b) The same as in panel (a), but normalized by the resistance 
of the junction in the normal state Rn- 



i.e. eI exc R N /A = (tt 2 /4) - 1. 



Discussion and conclusions. - In this work we have fo- 
cused our attention in the analysis of the I-V curves of 
short SNS graphene junctions, but the idea of combin- 
ing the transmission function of Eq. [21 with the single- 
channel point-contact results can be also used to study 
other transport properties of these junctions such as shot 
nois o 1 ^' 1 ? , Shapiro steps^S, photon-assisted tunneling 18 
or even full counting statisticsi2*2ii. For the special case 
of undoped graphene (V^ a t e = 0) the transport proper- 
ties are identical to those of short diffusive SNS systems, 
which have been already reported in the literaturei2iSS*2i. 
Moreover, following the recent work of Katsnelson^, a 
similar idea could be used to study superconducting con- 
tacts involving bilayer graphene. 

In summary, we have shown that the I-V characteris- 
tics of a short ballistic graphene SNS junction exhibits a 
very rich and gatable subharmonic gap structure. This 
very pronounced and sensitive structure predicted here 
can be straightforwardly compared with experimental re- 
sults in submicron scale junctions^ and it might consti- 
tute an ideal test quantity to confirm some of the pe- 
culiar consequences of the interplay of superconductivity 
and relativistic quantum dynamics. 

Acknowledgments.- We acknowledge discussions with 
Hubert B. Heersche, Pablo Jarillo-Herrero, Sebastian 
Bergeret and Elsa Prada. This work has been financed 
by the Spanish CYCIT (contract FIS2005-06255) and by 
the Helmholtz Gemeinschaft (contract VH-NG-029). 



1 K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. 
Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. Firsov, 
Science 306, 666 (2004). 

2 K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. 
I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and A. A. 
Firsov, Nature 438, 197 (2005). 

3 Y. Zhang, Y.-W. Tan, H. L. Stormer, and P. Kim, Nature 
438, 201 (2005). 

4 H.B. Heersche and P. Jarillo-Herrero, private communica- 
tion. 

5 C. W. J. Beenakker, Phys. Rev. Lett. 97, 067007 (2006). 

6 M. Titov and C. W. J. Beenakker, Phys. Rev. B 74, 041401 
(2006). 

7 P. Joyez, D. Vion, M. Gotz, M. H. Devoret, and D. Esteve, 
J. Supercond. 12, 757 (1999). 

8 T.M. Klapwijk, G.E. Blonder, and M. Tinkham, Physica 
B 109&110, 1657 (1982). 

9 E.N. Bratus, V.S. Shumeiko, and G. Wendin, Phys. Rev. 
Lett. 74, 2110 (1995). D. Averin and A. Bardas, Phys. 
Rev. Lett. 75, 1831 (1995). 

10 J.C. Cuevas, A. Martm-Rodero, and A. Levy Yeyati, Phys. 
Rev. B 54, 7366 (1996). 

11 E. Scheer, P. Joyez, D. Esteve, C. Urbina, and M.H. 
Devoret, Phys. Rev. Lett. 78, 3535 (1997). E. Scheer, 
N. Agrait, J.C. Cuevas, A. Levy Yeyati, B. Ludoph, A. 



Martm-Rodero, G. Rubio, J.M. van Ruitenbeek, and C. 
Urbina, Nature 394, 154 (1998). 

12 A. Bardas and D.V. Averin, Phys. Rev. B 56, R8518 
(1997). 

13 J. Tworzydlo, B. Trauzettel, M. Titov, A. Rycerz, and C. 
W. J. Beenakker, Phys. Rev. Lett. 96, 246802 (2006). 

14 Yu. V. Nazarov, Phys. Rev. Lett. 73, 134 (1994). 

15 S.N. Artemenko, A.F. Volkov, and A.V. Zaitsev, Zh. Eksp. 
Teor. Fiz. 76, 1816 (1979) [Sov. Phys. JETP bf 49, 924 
(1979)]. 

16 J.C. Cuevas, A. Martm-Rodero, and A. Levy Yeyati, Phys. 
Rev. Lett. 82, 4086 (1999). 

17 Y. Naveh and D.V. Averin, Phys. Rev. Lett. 82, 4090 
(1999). 

18 J.C. Cuevas, J. Heurich, A. Martm-Rodero, A. Levy Yey- 
ati, and G. Schon, Phys. Rev. Lett. 88, 157001 (2002). 

19 J.C. Cuevas and W. Belzig, Phys. Rev. Lett. 91, 187001 

(2003) . 

20 G. Johansson, P. Samuelsson, and A. Ingerman, Phys. Rev. 
Lett. 91, 187002 (2003). 

21 J.C. Cuevas and W. Belzig, Phys. Rev. B 70, 214512 

(2004) . 

22 M.I. Katsnelson, cond-mat /0606611 



